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Abstract: Very high energy (VHE, i.e. > 10 GeV) photons from Gamma-Ray Bursts (GRBs), as high as 90 GeV 
in rest frame energy, have been detected by the Fermi Large Area Telescope (LAT). This provides hope for a high 
statistics GRB detection by a ground-based instrument in the VHE domain. We here report on our estimates of the 
expected GRB detection rate and 7-ray rate (in case of detection) for the next generation ground-based Imaging Air 
Cherenkov Telescopes (lACTs): the Cherenkov Telescope Array (CTA). Moreover, we investigated the effect critical 
design parameters of CTA have on our predictions and we compare the performance of various observing strategies for 
the array pointing. Our estimations are based on a phenomenological model which uses temporal and spectral information 
from GRBs detected by Fermi-LAI as well as other instruments operating at lower energy. While detection of VHE 
emission from GRBs has eluded ground-based instruments thus far, our results suggest it should be within reach of 
CTA with a rate between 0.35 and 1.6 GRBs/year depending on the characteristics of the true GRB population and the 
performance that CTA will eventually achieve. 
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1 Introduction 

Detecting VHE emission from GRBs with ground-based 
telescopes would be a significant step forward not only 
for our comprehension of the underlying physics of these 
catastrophic events but also help us probe the evolving 
infrared-optical-ultraviolet extragalactic background light 
(EEL) at relatively high redshifts. Such VHE detection 
would most likely yield high photon statistics opening the 
door to detailed time-resolved spectroscopy in the VHE 
regime which would provide important clues on the emis- 
sion processes, particle distributions, and internal absorp- 
tion since current lack of statistics in this regime (only 
a handful of 7-ray photons were detected above 10 GeV 
by Fermi-LAI) prevent a clear distinction between differ- 
ent spectral models (power-law, broken power-law, smooth 
curvature, exponential cutoff...). Furthermore, by prob- 
ing the highest 7-ray energies emitted during prompt and 
afterglow emissions, VHE observations would set signif- 
icant constraints on the particle acceleration mechanisms, 
the possible presence of very-high energy cosmic-rays and 
the bulk Lorentz factor of the jet. Finally, by searching for 
signature of absorption through intervening IR-optical-UV 
light, important observational constraints on the EEL con- 
tent could be derived. 



GREs are a prime objective of current lACTs (MAGIC, 
VERITAS, HESS) but attempts of detection have proved 
fruitless so far. The next generation instrument CTA[jprob- 
ably holds the best chances at succeeding in this endeavor 
thanks in particular to its order of magnitude higher ef- 
fective area and few tens of GeVs energy threshold. A 
serendipitous observation of a GRE onset in the Field- 
of-View (FoV) of an lACT is extremely unlikely (~ 1% 
chance every year assuming a whole sky rate of ~ 600 
GREs/year and ~ 4° FoV). Therefore, Cherenkov tele- 
scopes rely on external GRE alerts in order to point their 
telescope arrays toward a GRE candidate for VHE emis- 
sion. We now make a little tour of GRE -dedicated mis- 
sion (to the best of our knowledge) that might be in orbit at 
the start of CTA operation currently scheduled for 2018, 
and would therefore be capable of providing such external 
alerts: 

• Swift Q: alert rate ~ 95 GREs/year with extremely 
good localization (< 10") 

• SVOM 11: alert rate - 70 - 90 GREs/year with 
similarly good localization; expected launch date 
- 2015 

1 . http://www.cta-observatory.org/ 
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• Fermi-GBM f3\: alert rate ^ 250 GRBs/year with 
poor localization (several degrees) 

Because the composition of the satellite fleet that will be 
in orbit at the time of CTA's first light remains quite uncer- 
tain, our approach is to provide independent estimates for 
the two vastly different type of GRB alerts that might be 
available: 1) Swift and/or SVOM, 2) Fermi-GBM. 

2 Modeling the GRB population at VHE 

'Make things as simple as possible, but no simpler' is the 
spirit we tried to follow to build our phenomenological 
model. The lack of knowledge about GRB VHE emis- 
sion naturally dominates the uncertainties in our predic- 
tions. We estimated the level of this uncertainty by consid- 
ering GRB populations with different spectral and temporal 
properties in the VHE domain, a range of behavior that can 
be reasonably assumed to encompass the properties of the 
true GRB population. 

2.1 Temporal and spectral model 

The VHE light curve is considered flat during the prompt 
emission to which a duration drawn from the BATSE Tgo 
distribution is assigned. Following the prompt emission, 
an extended VHE emission was modeled with a temporal 
decay t~^-^ similar to the decay measured on bright LAT 
GRBs il0 

Because the VHE spectrum is highly unknown and will 
strongly affect our predictions, we consider two different 
spectral models: 

• 'Bandex' model: a simple extrapolation of the Band 
function to VHE where we impose a maximum value 
of -2.0 to the high energy index. 

• 'Fixed' model: a power-law component is added on 
top of the Band function with an index -2.0 and nor- 
malization chosen to fix the energy flux ratio be- 
tween the LAT (lOOMeV - 300GeV) and BATSE 
(50 - SOOkeV) to 10%. 

Band function parameters are all drawn from the BATSE 
distributions. For Swift simulated bursts, a global fluence 
multiplier of 0.75 is applied to the BATSE fluence distribu- 
tion to provide the best fit between BATSE and SwiftfBAI 
fluences in the 15 — 150 keV range 

These two spectral scenarios are roughly consistent with 
the LAT detection rate of ~ 10 GRBs/year ('Bandex': 
- 10 GRBs/year; 'Fixed': - 20 GRBs/year). We also 
note that 'Bandex' and 'Fixed' spectral model match the 
behavior observed on the bright LAT GRBs: 080916C and 
090902B respectively. 

Let us point out that internal spectral curvature is not con- 
sidered mostly because of the large uncertainty this fea- 
ture creates. Indeed, depending on the level of curvature 



assumed, the VHE signal could either completely disap- 
pear or be almost unaffected. As a consequence of this 
caveat, the actual CTA observations could be substantially 
lower than our predictions if strong internal curvature be- 
low '^100 GeV are common in GRB spectrum. However, 
we point out that we consider very unlikely for the CTA 
GRB detection rate to be significantly higher than our pre- 
diction. In that sense, our estimates can be seen as upper 
limits for the actual CTA GRB detection rate. 

2.2 Redshift distribution 

Redshifts of our simulated bursts are drawn from the Swift 
population of ~ 170 GRBs with measured redshiftqj 
GBM (SVOM) is (will be) sensitive to a sHghtly different 
redshift distribution but we consider this difference to be a 
second order effect in our predictions. 

The EBL model used in our analysis is Gilmore 
& Somerville 2011 |5| approximately consistent with 
Gilmore et al. (2009), Franceschini et al. (2008), Finke et 
al. (2010), and Dominguez et al. (2010). We are currently 
investigating the effect that different EBL models would 
have on our estimates. 

3 CTA instrument response 

At the time of writing of this proceeding, the performance 
of the Cherenkov Telescope Array is still largely uncertain, 
in particular its low energy response which is crucial for 
GRB observations since extragalactic events suffer strong 
EBL absorption in the VHE domain. For this reason, we 
refrain from accurately modeling the CTA response but in- 
stead consider two CTA performances which we believe 
should encompass the final true CTA performance: 

• CTA 'baseline': 4 central Large Size Telescopes 
(LSTs) with energy threshold Eth = 25 GeV; 25 
Medium Size Telescopes (MSTs). 

• CTA 'optimistic': 4 LSTs with Eth = 10 GeV; 75 
MSTs. Background rate artificially lowered by a fac- 
tor of 3 to consider possible improvement in cosmic- 
ray rejection. 

Effective areas (see fig. [T]) and cosmic-ray rates for the 
LSTs and MSTs are derived through simple scaling of the 
well known VERITAS instrument performance and assum- 
ing a cosmic-ray spectrum ofE~'^-^. Small Size Telescopes 
(SSTs) are not considered since they are mostly sensitive 
to > 1 TeV showers. Zenith dependence of the energy 
threshold is assumed to be Eth = Eth,o x cos{Qzenith)~^ - 
As for the time delay it will take CTA to start observ- 
ing the field of a GRB from the time it was actually trig- 

2. temporal decay of t~^^ and t~^^ were also simulated with- 
out any significant change in the conclusions of our study. 

3 . http://heasarc.gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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Table 1 : ORB yearly detection rates for various lACT per- 
formance (MAGIC, CTA baseline, CTA optimistic), type 
of satellite alerts (Swift, GBM), spectral model ('bandex', 
'fixed'). Note that the GBM numbers are obtained when 
simulating the 'orbit' mode observing strategy (see section 
1431). 



Figure 1: Effective areas for MAGIC, VERITAS and the 
two CTA performances we are using in our study: baseline 
and optimistic. 

gerecj^ we assume 60 seconds for the LSTs and 100 sec- 
onds for the MSTs and show in section 4.2 the effect vary- 
ing this parameter has on our predictions. Finally our es- 
timates are based on one Cherenkov array only. In case 
the CTA collaboration is able to build a second array (one 
in each hemisphere) with similar characteristics, our pre- 
dicted rates would naturally increase by a factor of 2. 
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4 Results 

4.1 Detection and 7-ray rates 

For each of the GRB population models ('bandex' and 
'fixed') and each of the CTA performances ('baseline' and 
'optimistic'), we randomly simulated 10,000 GRBs in a 
solid angle 75° around zenith. We then computed the de- 
tection efficiency for these bursts, i.e. the probability of 
detecting a GRB when the Cherenkov array was able to 
observe the GRB field (i.e. for a GRB that occurs during 
a moonless night with good weather and at least 15° above 
the horizon). Assuming a Swift alert rate of 95 GRBs/year 
with a solar anti-bias factor of 1.4 |6| and a 10% duty cy- 
cle, table [TJprovides our predictions for the yearly GRB de- 
tection rate for various I ACT performances (MAGIC, CTA 
baseline, CTA optimistic), type of satellite alerts (Swift, 
GBM), and spectral model ('bandex', 'fixed'). In case such 
detection is established, fig. [2]shows it will almost certainly 
yield high photon statistics with a median number of 7-rays 
expected around 200-600 depending on the assumptions on 
the GRB model and the CTA performance. 

The redshift probability distribution of detected bursts is 
naturally shifted from the Swift distribution towards lower 
redshift due to the absorption by the EBL. The median red- 
shift of detected bursts is for example expected to be around 
redshift of 1.0. 



Figure 2: Number of 7-ray photons detected by CTA in 
case of detection 



4.2 Effect of critical CTA design parameters 

Figures |3] and [4] show how the predicted detection rates 
change with the two most critical parameters that we find: 
the LST energy threshold and time delay. Not too sur- 
prisingly, the energy threshold is crucial to catch the GRB 
emission which is softened by absorption through interac- 
tion with the EBL. Our study shows that a decrease in en- 
ergy threshold from ~ 25 GeV to ~ 10 GeV threshold 
would improve the GRB detection rate by a factor of ~ 2. 
The time delay for starting data taking has also a significant 
impact (fig. |4]). The CTA telescope slewing times that are 
aimed at for the design (LSTs: ~ 20 sec; MSTs: ~ 60 s) 
are already quite ambitious. A crucial aspect though will be 
to make sure that the other two factors contributing to the 
overall time delay (delay of the GCN alert to be received 
and time for the observers to authorize slewing) are as low 
as possible at the time CTA starts operation. 

We note that moonlight observations will probably be rou- 
tine for CTA therefore increasing the duty cycle from 
~ 10%to~ 13%. However this does not translate quite 



4. this time includes the time for the GCN alert to be sent out 
to the ground and received on site, for the observers to authorize 
slewing of the telescopes and for the actual slewing time of the 
telescopes to reach the target field. 
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Figure 3: Detection rate as a function of the energy thresh- 
old of the Large Size Telescopes. Our standard CTA 
baseline and optimistic performance assume 25 GeV and 
10 GeV respectively. 
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Figure 4: Detection rate as a function of the time delay for 
the Large Size Telescopes to reach the field of the GRB. 
Our standard CTA model assumes 60 seconds. 



linearly into a GRB detection rate increase since the energy 
threshold of such observations will be substantially higher. 

4.3 Optimizing CTA observing strategy for GBM 
alerts 

GBM alerts are a factor ^ 2.5 more frequent than Swift 
alerts. However the statistical and systematic uncertainties 
on the burst localization are on order of several degrees, 
and in most cases the error box will be substantially larger 
than the LST Field-of-View (FoV) targeted around 4° di- 
ameter. The simplest strategy for observing GBM bursts 
would be to point the CTA array at the best GBM local- 
ization and hope for the burst position to fall within the 
CTA FoV. However this results in extremely small detec- 
tion rates: ~ 0.05 GRB /year. A smarter strategy would 
be to implement a so-called 'orbit' mode where the GBM 
error box is being continuously scanned by the LST FoV. 
Predicted detection rates for such 'orbit' mode are signif- 
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Figure 5: Detection rate as a function of the field-of-view 
of the Large Size Telescopes. 

icantly more promising reaching a similar level as the de- 
tection rates from Swift alerts. Figure[5] shows the predicted 
detection rates for different assumptions on the LST FoV. 
We point out that for the 'orbit' mode, CTA will only be 
observing the field of the GRB a fraction of the time, there- 
fore only providing useful data on small chunks of the GRB 
light curve. A partial solution would be to achieve real-time 
detection and localization of the burst to subsequently point 
permanently on the GRB field. 

Figure[5] shows a very small improvement in detection rates 
when the LST FoV is varied which is due to the fact that 
the GBM error box is typically much larger. Therefore, 
a further increase in detection rate (as well as light curve 
coverage) of GBM bursts would necessitate a significant 
reduction of the GBM localization errors. Another observ- 
ing strategy has been proposed to cover the GBM error 
box: split the CTA array in multiple sub-arrays in order 
to achieve a uniform spatial and temporal coverage of the 
search region (which would come at the cost of lower sen- 
sitivity). However, this would only be possible with the 
MSTs since the LSTs are not numerous enough to cover 
the GBM error box. The energy threshold of MSTs being 
significantly higher, this observing mode would certainly 
lead to a much lower GRB detection rate than the 'orbit' 
mode which is therefore highly preferred by our study. 
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